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CHAPTER  1 


INTRODUCTION 


The  following  research  was  conducted  as  part  of  a  joint  project  for  the  University  of  Illinois 
and  the  Metallurgy  and  Quality  Assurance  Team  of  the  United  States  Army's  Construction 
Engineering  Research  Laboratory.  The  ultimate  objective  of  the  program  is  to  develop  an 
automated  welding  system  which  produces  consistently  good  welds  despite  variations  in  material 
parameters  and  other  disturbances.  The  scope  of  the  research  is  limited  to  consumable-electrode 
gas  metal  arc  welding  (GMAW)  since  this  is  one  of  the  most  frequently  employed  and 
economically  important  welding  processes. 

It  is  well  known  that  the  reliability  of  the  weld  is  strongly  correlated  to  the  microstructure 
and  gross  geometry  of  the  joint.  These  properties  are  determined  by  the  thermal  and  mechanical 
history  of  the  weld  puddle  and  the  rate  at  which  it  cools.  The  thermomechanical  dynamics  are 
driven  by  the  flow  of  heat  and  mass  from  the  torch  as  it  travels  along  the  weld  joint.  Arc  current  is 
one  of  the  key  process  variables  that  must  be  controlled  if  weld  quality  is  to  be  maintained.  This 
aspect  of  weld  quality  control  is  the  focus  of  this  work. 

A  gearmotor  drives  a  set  of  pinch  rollers  which  feed  wire  down  through  the  torch  head  for 
consumption  in  the  welding  process.  The  arc  current  is  determined  by  the  rate  of  wire  feed.  Arc 
current  control  is  accomplished  through  a  two-step  design.  First,  a  closed-loop  wire  feed  servo  is 
developed  which  can  track  a  reference  wire  feed  rate.  Second,  an  arc  current  controller  is  designed 
which  provides  the  reference  for  the  wire  feed  servo. 

Chapter  2  develops  the  wire  feed  mechanism  model  for  use  in  the  wire  feed  servo  controller 
design  of  Chapter  3.  Chapter  4  derives  the  arc  dynamics  model  which  relates  arc  current  to  wire 
feed  rate.  Chapter  5  makes  use  of  this  model  along  with  the  wire  feed  servo  dynamics  for  the 
closed-loop  arc  current  controller  design.  The  Appendix  provides  plots  derived  from  computer 
simulation  as  well  as  from  experimental  data  to  support  the  viability  of  the  design. 
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CHAPTER  2 


WIRE  FEED  MECHANISM  MODEL 


2.1.  Overview  of  the  Wire  Feed  System 


The  electro-mechanical  wire  feed  system  shown  in  Figure  2. 1  forms  the  plant  around 
which  the  closed-loop  wire  feed  servo  is  designed.  The  input  of  this  plant  is  taken  to  be  a  discrete¬ 
time  signal  u(k)  which  is  zero-order-held  ana  power-amplified  to  drive  a  dc  motor.  The  motor,  in 
turn,  drives  a  set  of  pinch  rollers  through  a  reduction  gearbox.  The  discrete-time  output  ex(k)  is 
defined  to  be  the  sampled  voltage  of  a  tachometer  on  the  motor.  This  signal,  properly  scaled, 
yields  a  sampling  of  the  wire  feed  rate  wp(t)  in  inches-per-minute. 


Torque  Disturbance 


Figure  2.1.  Wire  feed  mechanism  block  diagram. 

The  zero-order- hold  (ZOH)  digital  to  analog  converter  has  an  integer  input  u(k)  e  {-2048, 
...  ,  2048}  and  a  real  output  ec(t)  e  [-5,  5]  volts.  The  power  amplifier  is  modeled  as  a  fixed  x4 
gain  block.  This  simplification  of  the  power  amplifier  dynamics  is  valid  because  the  dominant  time 
constant  of  the  motor  is  two  orders  of  magnitude  larger  that  that  of  the  amplifier.  A  detailed 
analysis  of  the  dc  motor  dynamics  will  be  presented  in  Section  2.2.  The  system  employs  a  65.5: 1 
gearmotor  with  a  torque  efficiency  of  66%.  The  tachometer  conversion  constant  is  2  volts  per 
1000  rpm.  Noiseless  tach  measurements  are  assumed.  The  pinch  rollers  have  a  diameter  of 
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1.626  inches.  A  torque  disturbance  is  introduced  due  to  eccentricities  in  the  wire  roller  mechanism 
and  due  to  the  friction  of  the  wire  feed  path  as  shown  in  Figure  2.2. 


Figure  2.2.  Wire  path. 


2.2.  Direct  Current  Motor  Model 

The  dc  motor  of  Figure  2.3  is  modeled  as  a  series  connection  of  equivalent  armature 
resistance  and  inductance,  and  an  ideal  motor  (Kuo  (1980)). 


Figure  2.3.  Direct  current  motor  model. 


Signals  and  Constants: 


ea(t) 

armature  voltage 

Ra 

armature  resistance 

ia(t) 

armature  current 

La 

armature  inductance 

eb(t) 

back  emf  voltage 

Kb 

back  emf  constant 

®m(0 

shaft  velocity 

Jm 

rotor  inertia 

Tm(t) 

torque  developed  by  motor 

Ki 

torque  constant 

TL(t) 

load  torque 

Bm 

viscous  frictional  coefficient 

The  dc  motor's  dynamic  equations  are  now  developed.  Kirchoffs  Voltage  Law  around  the 
armature  loop  yields  the  electrical  constraint  equation 

ea(0  =  Ra  ia(t)  +  La^  +  eb(t) .  (2.1) 

Newton's  Second  Law  gives  the  mechanical  constraint  equation 

Tm(t)  =  TL(t)  +  BrnOWt)  +  (2.2) 

The  electro-mechanical  coupling  equations  are 

Tm(t)  =  Ki  ia(t) ,  and  eb(t)  =  Kb  com(t) .  (2.3) 

The  state  equations  for  the  motor  dynamics  are  now  obtained  by  substituting  equations  (2.3)  into 
(2.1)  and  (2.2). 


dia(t) 

dt 

dcouO) 

dt 


-  r1  ia(t)  -  ^  COm(t)  +  y— ea(t) , 


J^WmCO  +  f1^)  -  j~TL(t) 

J  m  J  m  J  m 


(2.4) 

(2.5) 


where  the  motor  constants 


Ra  =  2.55  x  101  Q 
La  =  2.04  x  10-3  H 


Kb=  2.31x10-2^ 

Jm  =  5.37  x  10'4  oz  in  s2 


Ki 

Bm 


3.27  x  101 


ozin 

~A~ 


1.00  x  10-3 


ozins 

rad 


are  obtained  from  the  manufacturer's  specification  sheets. 
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Analysis  of  the  system  time  constants  reveals  a  valid  simplification.  Make  the  following 
two  definitions: 

Electrical  time  constant  te  :  =  =  8.00xlO4s, 


Ra 

Mechanical  time  constant  xm  :=  Tp-  =  5.37  x  101  s. 

Time  constant  tm  exceeds  te  by  more  than  two  orders  of  magnitude.  This  suggests  that  it  is 
reasonable  to  neglect  the  eltctrical  dynamics.  Consider  (2.4)  multiplied  by  the  armature 
inductance  La : 


dia(t) 

dt 


—  "  Ra  ia(0  *  Rb  ®m(0  Jr  ®a(0  • 


(2.6) 


Since  La  is  small,  (2.6)  may  be  approximated  by 

0  —  -  Ra  ia(t)  -  Kb  com(t)  +  ea(t) . 

Solving  (2.7)  for  the  armature  current  ia(t)  gives 

ea(t)  -  Kb  com(t) 


(2.7) 


ia(0  - 


Ra 


(2.8) 


Substitution  of  (2.8)  into  (2.5)  yields  a  first-order  model  of  the  motor  dynamics: 

^  =  *  Pm  ®m(t)  +  j  Uj(t)  , 

where 


Bm  .  Ki  Kb 


Pm  :  ~  t  + 


Jm  Jm  Ra 


and 


UT(t)  :  =  ^  ea(t)  -  TL(t) 


(2.9) 

(2.10) 

(2.11) 


Taking  the  Laplace  transform  of  (2.9)  gives 

£2m(s) 


1  /  Jr 


UT(s)  S  +  pr 


The  resulting  transfer  function  description  of  the  dc  motor  is  shown  in  Figure  2.4. 
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TL(t) 


e  (t) 

a 

Ki 

u(t) 

1/Jm 

CO  (t) 
m 

Ra 

1  w 

S  +  Pm 

Figure  2.4.  Direct  current  motor  block  diagram. 


The  manufacturer's  specifications  for  the  wire  feed  motor  yield  the  following  numerical  values: 


{^=1.28^,  j-  =  1.86  x  103  — — j ,  pm  =  57.025 

Ka  V  ozin-s2  s 


Note  that  the  wire  feed  motor  acts  as  a  low-pass  filter  with  a  cutoff  frequency  of  9. 1  Hz. 


2.3.  Torque  Disturbance  Model 

The  wire  feed  gearmotor  turns  a  set  of  spring-loaded  pinch  rollers  which  pull  wire  from  a 
large  spool,  through  a  series  of  guides,  and  down  to  the  torch  head  (See  Figure  2.2).  The 
resulting  torque  load  on  the  motor  has  two  principal  components.  There  is  a  large  constant  load 
due  to  the  stiffness  of  the  wire  and  the  friction  of  the  guides,  rollers,  and  bearings.  In  addition, 
mechanical  imperfections  in  the  pinch  roller  assembly  cause  an  eccentricity  that  gives  rise  to  an 
approximately  sinusoidal  term. 

The  total  torque  load  on  the  gearmotor  shaft  is  modeled  as 

TG(t)  =  A  +  B  sin  (27t  fDt).  (2.12) 

This  torque  at  the  gearbox  is  a  product  of  the  torque  load  on  the  motor,  the  torque  efficiency  of  the 
gearbox,  and  the  gear  reduction  ratio.  Reflect  the  gearbox  torque  back  through  the  reduction  gears 
to  obtain  the  torque  load  on  the  motor  shaft: 

TL(t)  =  (0.66)1  (65.5)1  Tc(t)  =  0.0231  Tc(t).  (2.13) 
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The  frequency  of  the  sinusoidal  component  is  proportional  to  the  wire  feed  rate: 


fD  =  ~ 
n 


1  1  revolution  1  min 


1.626  in  60  s 


wp  =  3.26  x  10'3  wp  , 


(2.14) 


where  1.626  inches  is  the  diameter  of  each  pinch  roller  and  the  wire  feed  rate  wp  is  given  in  inches 


per  minute. 

Equations  (2.12)  and  (2.14)  reflect  a  key  assumption.  Wire  feed  rate  wp  is  taken  to  be 
constant  in  (2.14),  resulting  in  a  constant  disturbance  frequency  fo  in  (2.12).  However,  wp 
actually  fluctuates  in  response  to  computer  commands  during  control  of  the  weld  process.  This 
complication,  along  with  the  fact  that  Tl(0  can  also  affect  wp  through  the  motor  dynamics,  cause 
the  disturbance  to  be  not  truly  sinusoidal.  The  frequency  fo  actually  approximates  only  the 
principal  component  of  the  disturbance.  The  simplification  is  validated  by  the  fact  that  this 
principal  component  is  the  only  one  that  has  a  severe  effect  on  the  wire  feed  rate.  Figure  A.2  of  the 
Appendix  shows  a  typical  wire  feed  rate  signal  taken  from  experimental  open-loop  data.  A  Fourier 
transform  of  the  signal  shows  that  a  principal  component  is  indeed  present  at  near  the  predicted 
frequency. 

Extensive  open-loop  experiments  show  that  typical  values  for  constants  A  and  B  in  (2.12) 
are  141  and  15.0  oz-in  respectively.  Reasonable  wire  feed  rates  lie  in  the  200-400  IPM  range, 
resulting  in  disturbance  frequencies  of  0.65-1.3  Hz.  These  facts  yield  a  typical  torque  load  at  the 
motor  shaft  of 


Tl(0  *  3.27  +  0.35  sin  (2n  fot)  oz-in,  0.65  Hz  <  fo  <  1-3  Hz.  (2.15) 

Section  2.2  shows  that  the  wire  feed  motor  responds  to  frequencies  exceeding  9  Hz.  Since 
wire  feed  rate  is  a  principal  variable  in  the  control  of  weld  quality,  a  successful  wire  feed  servo 
controller  must  attenuate  this  disturbance  as  much  as  possible. 
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2.4.  Discretization  of  the  Wire  Feed  Mechanism  Model 

Since  wire  feed  rate  control  is  to  be  accomplished  with  a  digital  computer,  a  discrete-time 
model  of  the  wire  feed  mechanism  is  necessary.  Insertion  of  the  first-order  motor  model  of  Figure 
2.4  into  the  wire  feed  mechanism  block  diagram  of  Figure  2. 1  yields  Figure  2.5. 


TL(t) 


Figure  2.5.  Wire  feed  servo  block  diagram. 


The  reference  wire  feed  rate  WR(k)  is  given  in  inches  per  minute.  Figure  2.5  provides  for  a  two 
degree  of  freedom  controller  structure  that  features  both  feedforward  and  feadback  filters. 

The  scaling  factor  Kw  converts  motor  angular  velocity  com  (rad-s-2)  into  wire  feed  rate 
wp  (IPM).  Figures  2.1  and  2.2  yield 


_  1  7tT.626  in  1  revolution  60  s  _  „  _  IPM  s 

w  ~  65.5  1  revolution  jn  rad  1  m^n  ~  rad 


The  factor  Kt  converts  motor  angular  velocity  0^  into  the  tachometer  voltage  ej. 


1  2  V  1  revolution  60  s  0„n  V- 

~  5  1000  RPM  2 k  rad  1  min  “  X  U  rad 


Finally,  Kp  converts  the  ..  iputer  control  voltage  u  into  torque  (oz-in). 


KE 


4K, 

‘  Ra 


5.129 


oz-in 

V 
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Discretization  of  the  continuous-time  dc  motor  transfer  function  is  accomplished  through 
the  use  of  the  hold  equivalence  method  (Franklin  &  Powell  (1980)): 


G(z)  =  (1  -  z-1) 


. s  (s  +  Vjj 

The  Z-transform  may  be  evaluated  by  using  the  transform  pair 


1  /Jr 


(s(s 


+  a )/ 


z(l-e-aT) 
(z-  l)(z-eaT) 


T  :  =  sample  period 


This  yields 


Make  the  following  definition: 


where 


and 


1 

1  -  exp(-pmTs) 

(2.16) 

■Jm  Pm 

z  -  exp(-pmTs) 

G(z)  :  = 

Kd 

(2.17) 

z  -  Pd  ’ 

Kd  =  T  n-[  1  -  exp(-pmTs)  ] 

Jm  Pm 

PD  =  exp(-pmTs) . 


A  discrete-time  model  of  the  entire  wire  feed  servo  system  is  shown  in  Figure  2.6. 


TL(k) 


Figure  2.6.  Discretized  wire  feed  servo  block  diagram. 
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Analysis  and  controller  design  are  simplified  by  introducing  the  normalized  model  of 
Figure  2.7. 


A 

TL(k) 


Figure  2.7.  Discretized  and  normalized  wire  feed  servo  block  diagram. 


where 

WR(k)  =  ^  wR(k) , 
WF(k)  =  ^  wF(k)  , 


TL(k)  =  KDKTTL(k). 


Hff(z)  =  Ke  Kq  Kt  Hff(z)  , 


Hfb(z)  =  Hfb(z)  , 


The  controller  design  is  carried  out  on  the  normalized  system  for  simplicity.  Then  the  actual 
controller  filters  HFf(z)  and  Hfb(z)  are  found  from  the  normalized  ones. 
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CHAPTER  3 


WIRE  FEED  SERVO  CONTROLLER  DESIGN 

3.1.  Design  Objectives 

The  wire  feed  servo  is  a  closed-loop  system  within  the  arc  current  control  loop  as  shown  in 
Figure  3.1.  It  is  important  to  consider  the  system  as  a  whole  while  formulating  the  design 
objectives  for  each  subsystem.  This  nested  control-loop  structure  allows  the  two  controllers  to 
share  the  burden  of  the  total  design. 


The  wire  feed  servo  must,  of  course,  exhibit  zero  steady-state  error  for  a  constant  reference 
signal  wr(1c).  In  addition,  the  wire  feed  plant  is  subjected  to  a  disturbance  torque  of  the  form 

Tl(0  =  A  +  B  sin  (2k  foO  oz  in  ,  0.65  Hz  <  I'd  <  1.3  Hz.  (3.1) 

This  disturbance  must  be  attenuated  as  much  as  possible  since  it  falls  within  the  passband  of  the  arc 
dynamics. 
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High  frequency  arc  noise  occurs  which  should  be  attenuated  around  the  loop  so  that  it  does 
not  interfere  with  tracking  of  the  reference  signal  iR(k).  For  this  reason  the  wire  feed  servo's 
bandwidth  will  be  limited  to  4.0  Hz.  A  damping  ratio  £  of  0.707  i  will  yield  a  reasonable  step 
response. 

Chapter  1  ended  with  the  development  of  the  discretized  and  normalized  wire  feed  servo 
block  diagram  of  Figure  3.2. 


A 

TL(k) 


Figure  3.2.  Discretized  and  normalized  wire  feed  servo  block  diagram. 


The  controller  transfer  functions  Hff(z)  and  Hfb(z)  must  be  designed  to  accomplish  the  objectives 
stated  above. 

The  sampling  interval  Ts  is  selected  according  to  a  practical  rule  given  in  Astrom  & 
Wittenmark  (1984).  The  sampling  rate  should  be  six  to  ten  times  the  closed-loop  bandwidth. 

4  Hz  BW  =>  24  Hz  <  ^  <  40  Hz  =>  0.025  s  <  Ts  <  0.042  s  .  (3.2) 

A  sampling  interval  Ts  of  0.04  second  will  be  used  for  data  conversion. 

Numerical  values  for  the  discretized  plant  parameters  Kd  and  po  may  now  be  calculated: 

Kd  =  r— r— f  1  -  exp(-pmTs)  ]  =  29.3159, 

Jm  Fm 

and 

pD  =  exp(-pmTs)  =  0.102182. 
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The  design  objectives  relating  to  reference  tracking  characteristics  may  be  conveniently  met 
by  selecting  a  desired  reference-to-output  transfer  function.  Consider  the  continuous-time,  second- 


order  transfer  function 


H  M  =  mn2 

Cl(  '  WR(S)  s2  +  2I>„S  +  0)n2  ’ 


H(jto)  has  logio  magnitude  given  by 


/ co2 o  r  o) 

Hci(jco)  dB  =  -  20  log10  A  /  (  1  )2  +  (  ~^)2  . 


The  zero  steady-state  tracking  error  requirement  is  achieved  because  |  Hcl(0)  |  dB  is  zero 
for  any  choice  of  C,  and  The  damping  ratio  £  can  be  specified  directly  and  the  natural  frequency 
con  can  be  chosen  to  obtain  a  4  Hz  bandwidth.  Define  coc  to  be  the  cutoff  frequency.  We  must 


coc  =  2k  fc  =  2k  (4  Hz)  =  8k  radsJ  , 


3  dB  =  |  Hcl(jcoc)  |  dB  =  -  20  log10  A  /  (  1  -  ^)2  +  (  2^)2 

V  0)n2  C0n 


Simplification  of  (3.6)  yields 


2  =  (  1-^1)2  + 


Substitution  of  the  desired 


c°fV? 


results  in 


COn  =  C0c 


Thus  the  natural  frequency  is  equal  to  the  cutoff  frequency  for  this  choice  of  damping  ratio. 


The  desired  continuous-time  transfer  function  is  now  given  by 


Hd(s) 


_ 631.655 _ 

s2  +  35.5431s  +  631.655  ‘ 


(3.9) 


The  zero-order-hold  equivalent  discretization  of  H(s)  is  found  with  the  aid  of  Table  3.1  of 
Astrom  &  Wittenmark  (1984). 


it  /  n  .  _  WF(z)  _  biz  +  b2 
'  ’  WR(z)  z2  +  aiz  +  a2 


(3.10) 


The  coefficients  bi,  b2,  ai,  and  a2  are  found  as  follows: 


CO 


=  con V 1  -  c2  =  17.7715, 


b,  =  1  -  af^y  +  p|  =  0.307238, 
co 


a  =  exp(-  C«nTs)  =  0.491221  , 


b2  =  a2  +  a 


(S^Y-P 
1“ 


\ 


=  0.189566, 


p  =  cos(coTs)  =  0.757800  , 


ai  =  -  2ap  =  -  0.744495  , 


Y  =  sin(coTs)  =  0.652487  , 


a2  =  a2  =  0.241298  . 


Substitution  of  these  numerical  values  into  (3.10)  yields 


WF(z)  0.307238  z  +  0.189566 

Hcl(Z)  :  WR(z)  "  z2  -  0.744495  z  +  0.241298  ' 


(3.11) 


The  objectives  for  the  wire  feed  servo  controller  design  are  now  reduced  to  the  attenuation 
of  the  torque  disturbance  and  the  achievement  of  the  desired  reference-to-output  transfer  function 
of  (3.11).  Figure  A.l  of  the  Appendix  shows  magnitude  plots  of  the  open-loop  wirefeed  plant 
(plot  1)  and  the  desired  reference-to-output  transfer  function  (plot  2). 


3.2.  Proportional-Sum  Controller  Design 


A  simple  proportional-plus-sum  digital  control  strategy  is  employed  to  achieve  the  design 
objectives  as  closely  as  possible. 

Hff(z)  =  KP  +  =  K?--z  (Ki  ~  -  Hfb(z)  =  1  .  (3.12) 

Insertion  of  this  controller  into  the  block  diagram  of  Figure  3.2  yields  Figure  3.3. 

A 

TL(k) 


Figure  3.3.  PS  controlled  wire  feed  servo  block  diagram. 


The  accumulator  action  in  the  feed-forward  controller  gives  zero  steady-state  tracking  error 
for  constant  reference  inputs  as  well  as  rejection  of  constant  disturbances.  However,  since  only 
the  two  design  parameters  Kp  and  Ks  exist,  an  exact  matching  of  the  desired  reference-to-output 
transfer  function  of  (3.11)  is  impossible.  The  desired  characteristic  equation  can  be  obtained  but 
then  no  degrees  of  freedom  will  remain  with  which  to  place  the  closed-loop  zero. 

The  reference-to-output  transfer  function  is  given  by 


Hd(z) 


Wf(z) 

Wr(z) 


_ KPz  +  (Ks  -  KP) _ 

(z  -  0.102182)(z  -  1)  +  KPz  +  (Ks  -  KP)  ’ 


(3.13) 


Define  the  characteristic  equation  (denominator  of  the  closed-loop  transfer  function)  to  be  A(z). 


A(z)  =  Z2  +  (KP-  1.102182)  z  +  (KS-KP  +  0.102182).  (3.14) 
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The  desired  characteristic  equation  of  (3.1 1)  is 


A(z)  =  z2  -  0.744495  z  +  0.241298  . 


(3.15) 


By  setting  (3.14)  equal  to  (3.15)  and  equating  coefficients  of  like  powers  of  z,  one  obtains 


KP  =  -0.744495  +  1.102182  =  0.357687  (3.16) 

Ks  =  0.241298  -0.102182 +  KP  =  0.496803. 


The  resulting  proportional- sum  wirefeed  controller  is 


Hff(z) 


0.357687  z  +  0.139116 
z  -  1 


(3.17) 


Hfb(z)  =  1. 

Plot  3  of  Figure  A.  1  in  the  Appendix  shows  the  magnitude  of  the  closed-loop  disturbance-to- 
output  transfer  function  resulting  from  use  of  the  PS  controller  of  (3.17).  The  upper  half  of  Figure 
A. 3  shows  a  wirefeed  experiment  utilizing  this  controller.  Comparison  with  Figure  A. 2  shows  a 
reasonable  amount  of  low  frequency  disturbance  attenuation  over  the  open-loop  experimental 
results. 


3.3.  Polynomial  Controller  Design 

The  proportional-sum  controller's  simple  structure  does  not  utilize  the  feedback  filter  and 
this  method  only  gives  the  designer  two  parameters  to  tune.  In  order  to  better  accomplish  the 
design  objectives,  the  more  sophisticated  controller  of  Figure  3.4  will  be  used. 

The  design  methodology  presented  in  this  section  treats  the  various  transfer  functions  as 
rational  polynomials  in  z.  This  polynomial  method  provides  a  step-by-step  way  to  construct  a 
higher-order  controller.  See  Chapter  10  of  Astrom  &  Wittenmark  (1984).  This  design  will  exactly 
match  the  desired  reference-to-output  transfer  function.  It  will  also  improve  low  frequency 
disturbance  rejection  by  placing  a  1  Hz  zero  in  the  disturbance-to-output  transfer  function. 
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A 


TL(k) 


Figure  3.4.  Polynomial  representation  of  wire  feed  servo. 


Throughout  this  section,  H(z)  will  represent  a  rational  transfer  function  in  z,  while  all  other 
capital  letters  will  represent  polynomials  in  z.  A  polynomial  in  z  is  of  the  form 


P(z)  =  anzn  +  a^z1*-1  +  ...  +  aiz  +  ao  , 


where  the  degree  of  P(z)  is  n,  and  P(z)  is  said  to  be  monic  if  an  is  unity.  With  this  notation  in 
mind,  the  following  definitions  are  presented: 


Plant: 

H(z) 

ii 

>1“ 

ii 

1 

z  -  0.102182  ’ 

(3.18) 

Model: 

Hm(z) 

Bm(z) 

0.307238  z  +  0.189566 

(3.19) 

Am(z) 

z2  .  0.744495  z  +  0.241298  ’ 

Controller: 

Hff(z) 

,  T(z) 

’  ‘  R(z)  ’ 

M  .  S(z) 

Hfb(z)  .  -  . 

(3.20) 

The  rational  transfer  functions  of  the  controller  must  at  least  be  proper  in  order  for  the  controller  to 
be  causal  and,  therefore,  implemen table.  Thus  the  polynomials  R,  S,  and  T  must  be  such  that 

deg(R)  >  deg(T)  >  deg(S) .  (3.21) 

Furthermore,  restrict  polynomials  A,  Am,  and  R  to  be  monic  and  assume  no  pole-zero  cancellation 
in  either  of  the  transfer  functions  H  or  Hm. 
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The  plant  and  model  polynomials  can  now  be  identified: 


A(z)  =  z-  0.102182,  B(z)  =  1.  (3.22) 

Am(z)  =  z2  -  0.744495  z  +  0.241298  ,  Bm(z)  =  0.307238  z  +  0.189566  .  (3.23) 


Examination  of  Figure  3.4  yields  the  following  transfer  functions: 

WF(z)  =  HffH  Wf(z)  =  H 

WR(z)  1  +  HffHfbH  ’  Tl(z)  1  +  HffHFbH 


Substitution  of  (3.18)  and  (3.20)  into  (3.24)  gives 

WF(z)  _  BT  WF(z)  _  BR 

WR(z)  AR  +  BS  ’  Tl(z)  AR  +  BS  ' 


(3.25) 


Zero  steady-state  response  to  constant  disturbances  as  well  as  to  1  Hz  disturbances  is 
required.  The  disturbance-to-output  transfer  function  of  (3.25)  must  have  zeros  at  these 
frequencies.  A  polynomial  of  the  form 

P(z)  =  (z  -  r  ei0)(z  -  r  e‘Je)  =  z2  -  2r  cos0  z  +  r2  (3.26) 


has  complex  conjugate  zeros  at 

z  =  r  exp(±  j9) . 

With  the  radius  r  set  to  one,  P(z)  has  its  zeros  on  the  unit  circle  at  angles  ±0.  Thus  P(z)  will  zero  a 
sinusoid  of  frequency  0.  A  1  Hz  disturbance  has  a  discrete-time  normalized  frequency  given  by 

0  =  2?t  fD  Ts  =  2ti(1  Hz)(.04  s)  =  0.2513  rad.  (3.27) 

Consider  a  controller  polynomial  R(z)  of  the  form 

R(z)  =  (z  -  l)[z2  -  2  cos(0.2513)  z  +  1]  .  (3.28) 
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From  (3.25)  it  can  be  seen  that  a  zero  of  R(z)  will  also  be  a  zero  of  the  disturbance-to-output 
transfer  function.  The  (z  -  1)  term  provides  a  zero  at  dc,  while  the  quadratic  term  provides  zeros  at 
1  Hz. 

The  reference-to-output  transfer  function  is  required  to  be  Hm(z).  A  comparison  of  (3.19) 
with  (3.25)  yields  the  restriction 


BT  _  Bm_ 
AR  +  BS  Am 


(3.29) 


Since  R(z)  has  been  specified  as  (3.28)  and  since  A,  B,  Am,  and  Bm  are  also  fixed,  only  S  and  T 
remain  free  to  satisfy  (3.29).  Recalling  that  B  is  unity  for  the  normalized  plant,  it  may  be 
eliminated  in  order  to  write 


T(z)  =  Bm(z)  =  0.307238  z  +  0.189566. 


(3.30) 


S(z)  =  Am(z)  -  A(z)R(z) 

=  (z2  -  0.744495  z  +  0.241298)  -  (z  -  0.102182)(z  -  l)(z2  -  1.93718  z  +  1)  .  (3.31) 

Unfortunately,  this  solution  results  in  a  noncausal  controller  since  restriction  (3.21)  is 
violated.  Notice  that  the  degrees  of  R,  S,  and  T  are  three,  four,  and  one,  respectively.  This 
problem  is  elegantly  solved  by  the  introduction  of  an  observer  polynomial  A0(z).  Requirement 
(3.29)  is  unaltered  by  the  modification 


BT  _  BmA0 
AR  +  BS  “  AmAo  ’ 


(3.32) 


for  any  polynomial  A0.  However,  we  now  have  the  relations 


T(z)  =  Bm(z)A0(z)  =  (0.307238  z  +  0.189566)  A0  . 


(3.33) 


S(z)  =  Am(z)A0(z)  -  A(z)R(z) 

=  (z2  -  0.744495  z  + 0.241298)  A0  -  (z  -  0.102182)(z  -  l)(z2  -  1.93718  z  +  1)  .  (3.34) 
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Careful  examination  of  (3.33)  and  (3.34)  reveals  that  any  monic  A0  of  degree  fwo  will  yield 


deg(R)  =  deg(S)  =  deg(T)  =  3,  (3.35) 

resulting  in  proper  third-order  feedforward  and  feedback  controllers.  The  observer  polynomial 
derives  its  name  from  the  fact  that  it  is  cancelled  out  of  the  reference-to-output  transfer  function  just 
as  the  dynamics  of  a  state-space  observer  would  be.  It  does,  however,  appear  in  the  denominator 
of  the  closed-loop  disturbance-to-output  transfer  function. 

WF(z)  _  BT  _  BmA0  WF(z)  _  BR  _  R 

WR(z)  AR  +  BS  AmAo  ’  Tl(z)  AR  +  BS  AmAo  ' 


The  poles  of  this  quadratic  can  be  used  to  further  improve  the  disturbance  rejection  without 
modifying  the  reference  tracking  characteristics.  Their  location  was  chosen  by  an  iterative  process 
of  simulation.  It  was  found  that  a  frequency  of  1  Hz  and  a  radius  of  0.6  yield  the  most  desirable 
transfer  function.  Equation  (3.26)  gives 

A0(z)  =  z2  -  2(0.6)cos(0.087t)  z  +  (0.6)2  =  z2  -  1.16230  z  +  0.36  .  (3.37) 


Substitution  of  (3.37)  into  (3.33)  and  (3.34)  completes  the  design. 


R(z)  =  1.000000  z3  -  2.937166  z2  +  2.937166  z  -  1.000000, 

S(z)  =  1.132553  z3  -  1.770668  z2  +  0.751647  z  -  0.015315, 

T(z)  =  0.307238  z3  -  0.167536  z2  -  0.109727  z  +  0.068244. 


The  resulting  polynomial  method  wirefeed  controller  design  is  given  by 


HFF(z) 


0.307238  z3  -  0.167536  z2  -  0.109727  z  +  0.068244 
z3  -  2.937166  z2  +  2.937166  z  -  1.000000 


HFb(z) 


3.686248  z3  -  5.763188  z2  +  2,446468  z  -  0.049847 
z3  -  0.545299  z2  -0.357140  z  +  0.222120 


(3.38) 


(3.39) 
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Plot  4  of  Figure  A.l  in  the  Appendix  shows  the  magnitude  of  the  disturbance-to-output  transfer 
function  for  this  wirefeed  controller.  The  lower  half  of  Figure  A. 3  shows  the  results  of  a  wirefeed 
experiment  using  the  controller  of  (3.39).  A  comparison  of  this  design  with  the  proportional-sum 
design  of  Section  3.2  shows  that  the  polynomial  controller  improves  low  frequency  disturbance 
rejection  while  worsening  performance  at  higher  frequencies.  This  result  can  be  predicted  by 
examining  the  frequency  responses  of  Figure  A.l. 

3.4.  Two  Degree  of  Freedom  Controller  Design 

By  forcing  a  zero  at  1  Hz,  the  design  of  Section  3.3  accentuates  higher  frequency 
components  of  the  torque  disturbance.  A  controller  that  is  capable  of  rejecting  the  strong  low 
frequency  components  without  amplifying  the  higher  frequency  ones  is  very  desirable.  To  this 
end,  temporarily  ignore  the  reference  input  wr(1c)  and  treat  the  wire  feed  servo  as  the  simple  SISO 
regulator  of  Figure  3.5. 


A 

TL(k) 


plant 


i 

-1 

Z'PD 

controller 


w'pOO 


Figure  3.5.  Wire  feed  regulator. 

Let  the  controller  G(z)  have  a  simple  proportional-sum  structure. 

G(z)  =  KP  +  =  Kpz-MKs-Kp)  (3.40) 

The  accumulator  action  will  guarantee  zero  steady-state  error  for  constant  torque  disturbances  just 
as  it  did  in  Section  3.2. 
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The  closed-loop  disturbance  transfer  function  is  given  by 


TT  ,_N  Wp(z)  Z  -  1  A 

L(Z)  "  Tl(z)"  “  (z  -  pD)(z  -  1)  +  KPz  +  (Ks  -  KP) 

_ z  -  1 _ 

z2  +  (Kp  -  po  -  l)z  +  (Ks  -  Kp  +  po) 

The  PS  controller  allows  arbitrary  placement  of  the  closed-loop  poles.  The  objective  is  to 
minimize  the  regulator's  response  to  disturbance  inputs  all  the  way  from  dc  out  to  the  Nyquist 
frequency.  The  effect  of  pole  location  on  the  magnitude  of  the  frequency  response  will  now  be 
considered  with  the  aid  of  Figure  3.6. 


Re(z) 


Figure  3.6.  Z-plane  pole-zero  plot  of  Hl(z). 


Let  d,  be  the  absolute  distance  from  the  ith  pole  or  zero  to  the  point  on  the  unit  circle  that 
corresponds  to  the  frequency  of  interest.  Let  the  ith  pole  or  zero  be  located  at  the  point 

z  =  x  +  jy  . 

This  distance  may  be  expressed  as 


d,  =  I  eJ9  -  x  -  jy 


(3.42) 
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Now  let  the  poles  of  Hl(z)  be  located  at  a  ±  jb.  We  have 


and 


Hl(z) 


_ z  -  1 _ 

(z  -  a  -  jb)(z  -  a  +  jb)  ’ 


(3.43) 


HL(eJ0) 


ei0 

-  1 

eJe  -  a  - 

jb 

eJ0 

-  a  +  jb 

d, 

dpi  dp2 


(3.44) 


Examination  of  (3.44)  in  conjunction  with  Figure  3.6  gives  the  designer  considerable 
intuition.  The  zero  at  DC  causes  (3.44)  to  increase  for  higher  frequencies.  Extending  this  idea,  we 
now  see  why  an  additional  zero  at  1  Hz  causes  the  accentuation  of  high  frequency  disturbances. 
Conversely,  the  poles  decrease  (3.44)  as  they  are  placed  farther  away  from  the  frequency  of 
interest.  Since  the  objective  is  to  push  the  disturbance  down  as  much  as  possible  across  the  entire 
frequency  range,  it  is  readily  seen  that  placing  both  poles  at  the  origin  is  the  best  solution.  A 
controller  design  that  places  all  of  the  closed-loop  poles  at  the  origin  is  called  a  deadbeat  response 
design  (Kuo  (1980)). 

Thus  we  have  a  desired  closed-loop  disturbance-to-output  transfer  function: 


HlM  = 


_ z  -  1 _ 

Z2  +  (Kp  -  pD  -  l)z  +  (Ks  -  Kp+  pD) 


(3.45) 


Matching  coefficients  of  like  powers  of  z  in  the  denominators  yields 

Kp  -  po  -  1  =  0 ,  Kj  -  Kp  +  po  =  0 .  (3.46) 


Recalling  that  po  is  0. 102182  and  solving  for  the  controller  gains  give 

KP  =  1.1022,  Ks  =  1.0000, 

and 


G(z) 


KPz  +  (Ks  -  Kp) 
z  -  1 


1.1022z  -  0.1022 
z  -  1 


(3.47) 

(3.48) 


Now  that  the  desired  torque  disturbance  controller  G(z)  has  been  designed,  we  can  return 
to  the  full  multi-input  problem. 
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A 

TL(k) 


Figure  3.7.  Full  wire  feed  servo  block  diagram. 

A  comparison  of  Figures  3.5  and  3.7  results  in  the  relation 

G(z)  =  Hff(z)  Hfb(z)  .  (3.49) 


Also  from  Figure  3.7, 


Hr(z) 


and 


Hl(z) 


WF(z) 

HppH 

WR(z) 

1  +  HffHfbH  ’ 

WF(z) 

H 

Tl(z) 

1  +  HffHFbH 

A  second  relation  is  obtained  from  (3.50)  and  (3.51): 

Hr(z)  =  H^z)  Hl(z)  . 


(3.50) 

(3.51) 


(3.52) 


The  feedforward  and  feedback  controllers  can  now  be  found  from  (3.49)  and  (3.52). 

Hff(z)  =  Hl_1(z)  Hr(z)  ,  Hfb(z)  =  Hpp'Kz)  G(z) .  (3.53) 


Recall  the  desired  reference-to-output  and  disturbance-to-output  transfer  functions: 


Reference: 


Hr(z)  = 


0.307238  z  +  0.189566 
z2  .  0.744495  z  +  0.241298 


Disturbance:  Hl(z)  = 


z  -  1 


(3.54) 

(3.55) 
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Substitution  of  (3.54)  and  (3.55)  into  (3.53)  obtains  the  controller  transfer  functions.  The 
resulting  two-degree-of-freedom  wirefeed  controller  is  given  by 


Hff(z) 


0.307238  z3  +  0,189566  z2 _ 

z3  I  1.744495  z2  +  0.985793  z  -  0.241298 


(3.56) 


Hfb(z) 


3.587452  z3  -  3.003382  z2  +  1.1 13222  z  -  0.080242 
z3  +  0.617001  z2 


Plot  5  of  Figure  A.l  in  the  Appendix  shows  the  magnitude  of  the  disturbance-to-output  transfer 
function  for  this  wirefeed  controller.  Figure  A.4  shows  the  results  of  a  wirefeed  experiment  using 
the  controller  of  (3.56).  A  comparison  of  the  open-loop  and  the  three  conn-oilers  presented  in  this 
chapter  is  given  in  Figure  A. 5.  The  two-degree-of-freedom  controller  yields  superior  disturbance 
rejection  at  all  frequencies. 
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CHAPTER  4 


ARC  DYNAMICS  MODEL 


4.1.  Overview  of  the  Gas  Metal  Arc  Welding  Torch 

The  wire  feed  servo  rotates  a  set  of  pinch  rollers,  which  force  the  wire  into  the  torch  head 
and  through  the  contact  tube  whereupon  the  wire  is  consumed  by  the  welding  process,  as 
illustrated  in  Figure  4.1.  The  wire  feed  servo  is  in  itself  a  feedback  controlled  system  which  is 
capable  of  delivering  wire  to  the  weld  process  at  a  controlled  wire  feed  rate.  For  the  purposes  of 
this  derivation,  the  wire  feed  rate  wp(t)  is  considered  to  be  the  input.  The  arc  dynamics  form  the 
plant  and  the  arc  current  ia(t)  is  taken  to  be  the  output. 


Figure  4.1.  GMAW  torch. 
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Assumption  4.1. 


H  =  la(t)  +  ls(t)  =  constant  (4.1) 

where  H  is  the  contact  tube  height,  la(t)  is  the  arc  length,  and  ls(t)  is  the  stickout  length.  This 
assumption  is  motivated  by  careful  examination  of  the  side  profile  of  the  weld  bead  when  the  arc  is 
abrupdy  extinguished.  Figure  4.2  shows  how  the  weld  material  piles  up  behind  the  arc  as  opposed 
to  under  it. 


Figure  4.2.  Side  profile  of  weld  bead. 
Figure  4.3  shows  the  electrical  model  of  the  arc  dynamics. 

RS  Ls 


Figure  4.3.  Electrical  model. 
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Signals  and 

Constants: 

va(t) 

arc  voltage 

vs 

open  circuit  voltage  of  arc  welder  power  supply 

ia(0 

arc  current 

Rs 

Thevenin  resistance  of  arc  welder  power  supply 

vcr(t) 

contact  tube  voltage 

Ls 

plus  cabling  resistance 

inductance  of  arc  welder  power  supply 

Ris(t)  resistance  of  stickout  length  (time  varying) 

Assumption  4.2. 

The  electric  arc  is  modeled  as  shown  in  Figure  4.4,  where  V0  is  the  sum  of  the  cathode  and 
anode  voltage  drops,  E0  is  the  electric  field  intensity  of  the  arc  plasma,  and  la(t)  is  the  length  of  the 
electric  arc. 


i(0 

a 


Figure  4.4.  Arc  model. 

There  are  two  modes  associated  with  the  arc: 

Te  =  electrical  mode;  the  arc  can  be  modelled  as  a  time  varying  conductance. 

Tt  =  thermodynamic  mode;  consumption  of  the  wire. 

Time  constant  Xj  exceeds  Te  by  approximately  two  orders  of  magnitude.  A  static  model  is 
therefore  used  for  the  arc  electrical  circuit.  However,  the  dynamics  of  wire  consumption  will  not 
be  ignored. 


28 


4.2.  Development  of  the  Dynamic  Equations 


Kirchhoff  s  Voltage  Law  applied  to  Figure  4.3  yields  the  electrical  constraint  equation, 

Vs  =  Rsia(t)  +  Ls^  +  Ris(t)  ia(t)  +  V0  +  E0la(t).  (4.2) 


The  stickout  length  resistance  R^t)  can  be  written  as  a  function  of  the  arc  length: 

Ris(t)  =  j  ls(0  =  ^  [  H  -  la(t)  ] 


(4.3) 


where  p  is  the  resistivity  of  the  weld  wire  (  Q  •  in  )  and  A  is  the  cross-sectional  area  of  the  weld 
wire.  Substitution  of  (4.3)  into  (4.2)  yields  the  electrical  constraint  equation  in  final  form: 


dia(t) 

dt 


-£-[Rs  +  £-H]ia(t)  +  ^iaWWt)  -  ^la(t) 


V0 


(4.4) 


Define  wF(t)  to  be  the  feed  rate  of  weld  wire  and  define  wc(t)  to  be  the  consumption  rate  of 
weld  wire  in  inches  per  minute.  Since  la(t)  is  the  length  of  arc  in  inches,  the  following  relation 
must  hold: 


dla(t) 

dt 


1  min  r 
~  60  s  1 


wc(t)  -  wF(t)]. 


(4.5) 


Following  Waszink  &  Van  Den  Heuvel  (1982), 

wc(t)  =  [  ki  ia(t)  +  k2  ia2(t)  ls(t)  ] ,  (4.6) 

where  the  second  term  is  due  to  joule  heating. 

Assumption  4.1  and  Equations  (4.5)  and  (4.6)  yield  the  mechanical  constraint  equation: 

“fp  =  k!  ia(t)  +  k2  H  ia2(t)  -  k2  ia2(t)  la(t)  -  ^cop(t) .  (4.7) 


The  arc  current  ia(t)  and  the  contact  tube  voltage  vct(0  are  measured  quantities,  while  the 
arc  length  la(t)  is  not  measured.  It  is  for  this  reason  that  an  additional  relation  is  of  interest. 
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Kirchhoff  s  Voltage  Law  applied  to  an  inner  loop  of  Figure  4.3  gives 


vct(0  =  Ris(0  ia(0  +V0  +  E0la(t) 


Substitution  of  (4.3)  into  (4.8)  allows  la(t)  to  be  measured  indirectly: 


vct(0  =-^ia(t)  -  ^ia(t)la(t)  +  E0la(t)  +  V0 


Summary  of  the  arc  dynamics  model: 


^a(t)  _  1  r  p  j.  £u  1  i  fA  x  P  :  i  i  /f\  .  Vs  ~  Vo 

^  _  "  L  A H  ^a(t)  *aw  ta(t)  -  ta(t)  + 


^  =  ki  ia(t)  +  k2Hia2(t)  -  k2ia2(t)la(t)  -  ^wp(t) 


(4.10) 


vcr(t)  =-^ia(t)  -  ^-ia(t)la(t)  +  E0la(t)  +  V0 


where 


ia(t)  =  arc  current, 
la(t)  =  arc  length. 


vcr(t)  =  contact  tube  voltage, 


wp(t)  =  wire  feed  rate. 


Vs  =  open  circuit  voltage  of  power  supply, 

Rs  =  resistance  of  power  supply  and  cabling, 
Ls  =  inductance  of  power  supply, 
p  =  resistivity  of  weld  wire, 

A  =  cross-sectional  area  of  weld  wire, 

H  =  contact  tube  height, 

V0  =  sum  of  cathode  and  anode  voltage  drops, 
E0  =  electric  field  intensity  of  arc  plasma, 
kj,k2  =  wire  consumption  rate  constants. 
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4.3.  Calculation  of  the  Parameter  Values 

The  following  nominal  operating  point  (measured  experimentally)  is  used  as  a  basis: 
ia  =  340  A.  vct  =  31  V. 

la  =  0.25  in.  coF  =  300  IPM. 

The  power  supply  is  set  for  a  Vs  of  40  V,  and  the  manufacturer's  specifications  show  an  L$  of 
400  p.H.  The  0.0625-inch  diameter  weld  wire  has  a  cross-sectional  area  A  of  3.068xl0‘3  in2. 
The  contact  tube  height  H  is  set  at  1  in.  Wu  &  Richardson  (1984)  approximate  p  at  40  |i£2in. 
Table  A.l  of  that  paper  yields  an  E0  of  26  V-in*1. 

For  the  calculation  of  V0  and  Rs,  consider  the  electrical  model  in  steady-state  as  shown  in 
Figure  4.5. 


Figure  4.5.  Electrical  model  in  steady-state. 


Analysis  of  the  circuit  in  steady-state  yields  the  following  relations: 

RS  =  Vs^VcT> 

V0  =  Vcr  -  E0la  -  Ris  ia  =  Vcr  -  E0la  -  £(H  -  la)  ia .  (4.11) 

Substitution  of  the  numerical  values  gives  an  Rs  of  0.02647  £2  and  a  V0  of  21. 18  V. 
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Estimation  of  the  constants  ki  and  k2  requires  use  of  experimental  data.  Curve  fitting  with 
actual  welding  runs  yielded  the  approximate  relation 

wc  =  0.59398  ia  +  0.00047164  ia2 ,  (4.12) 

where  wc  is  in  inches  per  minute  and  ia  is  in  Amps.  Setting  (4.12)  equal  to  (4.6)  gives 

0.59398  ia  +  0.00047164  ia2  =  60  (  kj  ia  +  k2  ia2  1* ) .  (4.13) 

Substitution  of  the  operating  point  values  for  ia  and  ls  into  (4.13)  gives 

ki  =  9.8997xl0-3  in-A^-s-1  ,  k2  =  1.0481xl0-5  A-2-s-i  . 

These  values  of  ki  and  k2  yield  a  wire  consumption  rate  wc  of  260  IPM  for  an  arc  current  of 
340  A.  This  rate  is  within  10%  of  the  experimentally  measured  operating  point  value. 

The  arc  dynamics  model  is  now  complete  with  numerical  coefficients. 

Arc  dynamics  model: 

=  -  9.8770X101  ia(t)  +  3.2595X101  ia(t)  la(t)  -  6.5xlOMa(t)  +  4.705X104 

(4.14) 

^- =  9.8997x10-3  ia(t)  +  1.048 lxlO-5  ia2(t)  -  1.0481x10-5  ia2(t)  la(t)  -  1.667xl0-2  wF(t) 

vcr(0  =  1.3038xl0-2  ia(t)  -  1.3038x10-2  ia(t)  la(t)  +  2.50x101  la(t)  +  2.118x10! 
with  ia(t)  in  amperes,  la(t)  in  inches,  vcr(t)  in  volts,  and  wp(t)  in  inches  per  minute. 


4.4.  Linearization  of  the  Model 

A  linear  approximation  to  the  arc  dynamics  model  is  used  for  the  controller  design.  Define 
input  u(t)  to  be  the  wire  feed  rate  wp(t).  Define  state  xj(t)  to  be  the  arc  current  ia(t),  and  define 
state  x2(t)  to  be  the  arc  length  la(t)  in  order  to  write  the  nonlinear  state  equations  of  (4.15). 
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*i(t)  =  Cl  Xi(t)  +  c2  Xi(t)  x2(t)  +  c3  x2(t)  +  di  :  =  /i(xi,  x2)  +  dj 

x2(t)  =  c4xi(t)  +  c5xi2(t)  +  c6  Xi2(t)  x2(t)  +  b2u(t)  :=  /2(xx,  x2)  +  b2  u  (4.15) 

Now,  make  some  further  definitions: 


’xf 

x  :  =  , 

Lx2  J 

r/i(xi»x2) 

•"  L/2(X1.X2)J’ 

cr 

II 

\ - 1 

> 

d  :  = 

rdii. 

.0  j’ 

(4.16) 

aij  :  = 

d/i(xl,  x2) 

A  :  = 

an 

ai2l 

(4.17) 

3xj 

• 

a2i 

^22]' 

Using  the  above  vector  notation,  the  system  becomes 

x  =  /(x)  +  bu  +  d.  (4.18) 

Furthermore,  the  linearized  system  is  given  by 

8x  =  A  8x  +  b  8u  ,  (4.19) 

where  the  partial  derivatives  of  the  A  matrix  are  evaluated  at  a  nominal  operating  point  xnom  and  the 
variables  8x  and  8u  represent  reasonably  small  perturbations  from  this  operating  point 
(Brogan  (1982)).  The  partial  derivatives  of  (4.15)  are  taken  to  obtain 

an  =  ci  +  c2x2,  ax2  =  c2xi  +  c3 , 

a2i  =  c4  +  2  C5  xi  +  2  C6  xi  x2  ,  a22  =  C6  xi2.  (4.20) 

Now,  since  the  arc  current  ia(t)  is  the  variable  to  be  regulated  and  the  arc  length  la(t)  is  not 
really  of  interest,  define  the  system  output  y(t)  to  be  the  arc  current  state  variable  xx(t).  In  vector 
notation, 

Sy  =  c  8x ;  c  :  =  [  1  0].  (4.21) 

Finally,  the  single-input,  single-output  system  of  (4.19)  and  (4.21)  has  a  transfer  function 
H(s)  that  is  given  by 

H(s)  =  c(sl  -  A)-!  b  .  (4.22) 


If  the  characteristic  polynomial  of  A  is  denoted  by  A(s),  then  H(s)  can  be  written  as 

H(s)  =  ;  A(s)  =  s2  -  (an  +  a22)  s  +  (an  a22  -  a]2  a21)  .  (4.23) 

A(s) 


Consider  three  values  for  the  arc  current;  200  A,  340  A,  and  400  A.  These  values  span 
the  range  of  reasonable  arc  currents  for  welding.  Now  consider  the  nonlinear  dynamic  equations 
of  (4.14)  in  steady-state  (derivatives  equal  to  zero).  These  equations  can  be  used  to  find  steady- 
state  values  for  the  arc  length  and  the  wire  feed  rate  that  correspond  to  each  arc  current. 


0  =  -  9.8770X101  ia  +  3.2595xl0Ha  la  -  b.SxlO4^  +  4.705xl04 


(4.24) 

0  =  9.8977xl0-3  ia  +  1.0481xl0-5  ia2  -  1.0481xl0'5  ia2  la  -  1.667xlO“2wF 


Substitution  of  the  arc  currents 
ia  =  200  A. 
ia  =  340  A. 
ia  =  400  A. 


into  (4.24)  yields  the  following  three  operating  points: 
=>  la  =  0.4668  in,  wF  =  132 IPM, 


=>  la  =  0..2498  in,  wF  =  256  IPM, 

=>  la  =  0.1451  in,  wF  =  324  IPM. 


(4.25) 


Substitution  of  the  numerical  values  of  (4.14)  and  (4.25)  into  the  relations  of  (4.15),  (4.20),  and 
(4.23)  yields  the  following  three  linearized  arc  dynamics  transfer  functions: 


,  _  _ 974,68 _ 

H(S)  200A  s2  +  83.975  s  +  744.72  ’ 


(4.26) 


H(s)  340 A 


_ 898.63 _ 

s1  +  91.839  s  +  931.86  ’ 


(4.27) 


H(s)  400 A 


_ 866.04 _ 

s2  +  95.700  s  +  1044.6  ' 


(4.28) 


Figure  4.6  shows  how  the  poles  of  the  linearized  transfer  function  H(s)  vary  (in  rad/sec) 
with  respect  to  nominal  arc  current.  This  variation  is  smooth  and  moderate,  suggesting  that  a 
single  H(s)  can  be  used  for  the  design  of  the  arc  current  controller.  Choose  this  nominal  transfer 
function  to  be  the  340  A  H(s)  of  equation  (4.27). 
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Figure  4.6.  Pole  locations  versus  nominal  arc  current. 
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CHAPTER  5 


ARC  CURRENT  CONTROLLER  DESIGN 


Chapter  4  ended  with  the  development  of  a  nominal  continuous-time  transfer  function  that 
relates  arc  current  to  wire  feed  rate.  Although  the  arc  dynamics  are  nonlinear,  it  was  shown  that 
the  poles  of  the  linearized  system  do  not  shift  radically  when  the  nominal  operating  point  is  varied. 
It  is  therefore  reasonable  to  use  this  nominal  transfer  function  for  the  arc  current  controller  design. 

Chapter  2  posed  the  reference-tracking  objectives  for  the  wire  feed  servo  in  terms  of  a 
desired  reference-to-output  transfer  function.  This  transfer  function  was  matched  by  the  controller 
that  was  designed  in  Chapter  3. 


Figure  5.1.  Closed-loop  arc  current  system. 

Figure  5.1  shows  the  arc  current  system  with  ia(t)  and  ip Ck)  given  in  amperes,  and  wp(t) 
and  WR(k)  given  in  inches  per  minute.  The  closed-loop  system  is  required  to  have  a  bandwidth  of 
1  Hz.  Furthermore,  the  closed-loop  step  response  should  be  without  overshoot.  The  system  is  to 
be  mildly  overdamped.  Zero  steady-state  tracking  error  is  required  for  constant  reference  inputs. 
The  arc  current  controller  will  be  designed  in  a  continuous-time  setting  and  then  discretized  for 
computer  implementation. 
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The  design  objectives  are  not  difficult  but  the  plant  is  nonlinear  with  uncertain  parameters. 
For  this  reason,  a  simple  and  robust  proportional-integral  controller  will  be  used  for  the  design. 
Robust  stability  of  the  arc  current  is  of  primary  importance  for  reasons  of  safety  as  well  as  weld 
quality. 


Figure  5.2.  Continuous-time  arc  current  block  diagram. 


Figure  5.2  shows  the  arc  current  system  in  a  continuous-time  format.  The  output  ia(t)  of 
the  arc  current  plant  and  the  reference  input  iR(t)  are  given  in  amperes.  A  "Hall  effect”  sensor 
produces  a  voltage  that  is  proportional  to  the  arc  current.  A  constant  Kh  of  1  V  per  100  A  models 
this  sensor.  The  constant  Kw/Kt  converts  volts  to  inches  per  minute.  Its  value  is  found  in 
Chapter  2  to  be  194.95  IPM  per  volt.  Finally,  the  wire  feed  servo  has  its  input  and  output  given  in 
inches  per  minute. 

Chapters  2  and  4  give  the  closed-loop  wire  feed  servo  and  the  open-loop  arc  current 
transfer  functions. 


Closed-loop  wire  feed: 


Hw(s) 


_ 631.655 _ 

s2  +  35.5431s  +  631.655  ' 


(5.1) 


Open-loop  arc  current: 


Ha(s) 


_ 898.631 _ 

s2  +  91.8393s  +  931.862  ‘ 


(5.2) 


The  proportional-integral  arc  current  controller  is  of  the  form 


Hpi(s) 


K  Ki  =  KpL+Jm 

r  c  c 


(5.3) 
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The  controller  design  is  simplified  by  introducing  the  normalized  model  of  Figure  5.3. 


A  A 


Figure  5.3.  Normalized  arc  current  system. 


where 

A  A 

iR(t)  =  KniR(t),  ia(t)  =  KniaW, 

and 

Hpi(s)  =  KP  +  =  £0^  Hpi(s) .  (5.4) 

The  composite  plant  Hw(s)Ha(s)  is  fourth  order  with  relative  degree  four.  To  this  we 
append  a  first-order  controller,  creating  a  fifth-order  closed- loop  system.  A  purely  analytic  design 
of  Hpi(s)  to  achieve  a  desired  closed-loop  bandwidth  and  damping  would  be  difficult.  For  this 
reason,  the  design  was  accomplished  with  the  aid  of  iterative  computer  simulation.  Integrator  gain 
was  adjusted  to  achieve  the  desired  bandwidth,  and  then  the  proportional  gain  was  tuned  to  create  a 
mildly  overdamped  step  response.  The  resulting  gains  are 

KP  =  0.29,  Ki  =  4.20.  (5.5) 

Thus  the  continuous- time  proportional-integral  controller  is  given  by 

~  4  20 

Hpi(s)  =  0.29  +  .  (5.6) 
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A  simple  Euler  approximation  with  a  sampling  period  Ts  of  0.04  second  is  used  for  the 
discretization  of  the  PI  controller.  The  resulting  proportional-sum  controller  is 


Hps(z)  =  KP  + 


KiTs 


=  0.29  + 


0.168 
z  -  1 


0.29z  -  0,122 
z  -  1 


(5.7) 


Figure  A. 6  of  the  Appendix  shows  an  experimental  welding  run.  The  arc  current  controller 
was  stepped  across  nearly  the  entire  operating  range  to  show  the  robust  quality  of  the  design.  The 
upper  plot  shows  the  arc  current  and  the  lower  plot  shows  the  closed-loop  wirefeed  reference  wr(0 
(dashed  line)  and  the  wirefeed  servo  output  wp(t)  (solid  line). 
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CHAPTER  6 


CONCLUSIONS 


From  an  academic  point  of  view,  the  work  exhibits  some  diversity.  In  Chapter  2  a  wire 
feed  plant  is  developed  which  has  simple  linear  dynamics  but  is  corrupted  by  a  severe  torque 
disturbance.  This  model  is  discretized  and  Chapter  3  presents  three  discrete-time  controller  designs 
which  vary  considerably  in  performance,  complexity,  and  design  philosophy.  Chapter  4  develops 
a  nonlinear  arc  dynamics  model  which  is  linearized  for  the  controller  design.  Finally,  in  Chapter  5 
a  continuous-time  arc  current  controller  is  designed  and  then  discretized  for  computer 
implementation. 

The  plots  in  the  Appendix  give  testimony  to  the  success  of  the  final  design.  The  torque 
disturbance  in  the  wire  feed  servo  has  been  virtually  eliminated  while  the  servo  exhibits  desirable 
reference-to-output  characteristics.  The  fixed-parameter  arc  current  controller  is  well  behaved  at  all 
operating  points  and  is  capable  of  stepping  the  arc  current  the  full  100  A  of  its  useful  range. 


40 


APPENDIX 


SELECTED  PLOTS 


This  appendix  presents  experimental  data  as  well  as  computer  simulated  results.  These 
plots  are  gathered  together  to  facilitate  a  convenient  performance  comparison  of  the  various 
controller  designs.  All  Fourier  transforms  are  presented  on  a  decibel  scale  and  have  been 
normalized  by  dividing  by  the  average  value  of  the  time  domain  signal. 

Figure  A.l  gives  various  magnitude  plots  that  pertain  to  the  wire  feed  servo  controller 
design.  In  numbered  order  these  plots  are:  the  open-loop  plant,  the  desired  reference-to-output 
transfer  function,  and  the  closed-loop  disturbance-to-output  for  the  proportional-sum,  polynomial, 
and  two-degree-of-freedom  controller  designs. 

Figure  A.2  shows  open-loop  experimental  data.  Time  domain  plots  and  resulting  Fourier 
transforms  are  given  for  the  wire  feed  rate  and  for  the  resulting  arc  current.  Figures  A.3  through 
A. 5  show  experimental  results  for  the  three  wire  feed  servo  controller  designs.  Figure  A. 6  shows 
experimental  step  response  data  for  the  final  closed-loop  arc  current  controller. 


41 


MAGNITUDE  (PE) 


HIRE  FEED  RATE  <IPM> 


Figure  A.2.  Open-loop  experimental  data. 
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WIRE  FEED  RATE  (1PM) 


Figure  A. 5.  Open-loop,  PS,  polynomial,  and  TDF  wirefeed  FFT  comparison 
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